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Abstract 

We consider the initial value problem for the thermal-diffusive combustion systems 
of the form: u\ t t = A x u\ — uiu™, U2,t = dA x U2 + u\u 2 n -, x £ R n , n > 1, m > 1, d > 1, 
with bounded uniformly continuous nonnegative initial data. For such initial data, 
solutions can be simple traveling fronts or complicated domain walls. Due to the 
well-known thermal-diffusive instabilities when d, the Lewis number, is sufficiently 
away from one, front solutions are potentially chaotic. It is known in the literature 
that solutions are uniformly bounded in time in case d < 1 by a simple comparison 
argument. In case d > 1, no comparison principle seems to apply. Nevertheless, we 
prove the existence of global classical solutions and show that the L°° norm of v>2 can 
not grow faster than 0(loglogt) for any space dimension. Our main tools are local LP 
a-priori estimates and time dependent spatially decaying test functions. Our results 
also hold for the Arrhenius type reactions. 



1 Introduction 



In this paper, we are concerned with the existence of global classical solutions and large 
time asymptotic bounds of the thermal-diffusive combustion system: 

u ht = A x ui - mu™, 

u 2 , t = dA x u 2 + uiu™, (f.f) 

with nonnegative initial data (u 1: u 2 )\ t=0 = (o-i(x), a 2 (x)) G (C® u (R n )) 2 , the space of uni- 
formly bounded continuous functions on R n . Here x G R n , n, m are positive integers, 
d > 1 is the Lewis number; is the n dimensional Laplacian. System (1.1) describes the 
evolution of mass fraction of reactant A, u±, and that of the product B, u 2 , for the autocat- 
alytic chemical reaction of the form: A + mB — > (m+ 1)B with rate k m U\U™ k m a positive 
constant. In case m — 1, or 2, we refer to Billingham and Needham, [5], [6], for details. 
System (1.1) also describes the mass fraction, u±, and temperature, u 2 , of reactant A, of 
a one step irreversible reaction A — > B; especially when u™ is replaced by the Arrhenius 
reaction term exp{— E/u 2 }, E > being the activation energy. In this context, system 
(1.1) is the well-known thermal diffusive system, see Matkowsky and Sivashinsky [19]. 

One of the basic questions for (1.1) with L°° initial data is the existence of global 
solutions and the possible uniform in time bounds of u 2 . In case of the Arrhenius reaction, 
i.e. with Uiexp{— ^} replacing U\U™ in (1.1), there are many works on global solutions, 
see Avrin [2] , Larrouturou [15] for results in one space dimension, among others. Yet their 
bounds of the solutions grow linearly in time. It is still a conjecture whether u 2 is bounded 
uniformly in time, see Berestycki and Larrouturou [3], and Manley, Marion, and Temam 
[16]. 

On the other hand, system (1.1) on a bounded domain Q in R n with homogeneous 
boundary conditions has been thoroughly studied. The problem on existence and uniform 
bounds of solutions was first posed by R. Martin, and later solved partially by Alikakos 
[1], and completely by Masuda [18]. See also Haraux and Youkana [13], Hollis, Martin, 
and Pierre [14] for related approaches and extensions. System (1.1) on the line with spa- 
tially decaying data in L 1 n L°°(i? 1 ) has been investigated recently in Berlyand and Xin 
[4], Bricmont, Kupiainen, and Xin [9] for critical nonlinearity m — 2. Global classical so- 
lutions exist for any size initial data and converge to self-similar solutions with anomalous 
exponents [9]. 

System (1.1) with L°° data on R n is very different from either the one on the bounded 
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domains or the one with spatially decaying data. The system admits propagating front 
solutions, from simple traveling wave solutions to the complicated domain walls. When 
m — 1, 2, existence of traveling fronts are proved in [5]; in [6], formal asymptotic as well as 
computational studies are done for fronts generated from initial data u\ = positive constant, 
u 2 = bounded nonnegative function with compact support. In the Arrhenius reaction and 
high activation energy limit (E — > +oo), it is well-known that planar fronts are subject 
to (thermal-diffusive) instabilities when d is far enough away from one, and fronts become 
chaotic, see Clavin [10], Sivashinsky [20] , Terman [21] and references therein. In the 
long wave asymptotic limit, the perturbations of the planar fronts satisfy the celebrated 
Kuramoto-Sivashinsky equation [20]. For an interesting study on stable and unstable planar 
fronts away from the large E limit, see Bonnet et al. [7]. In spite of the front instabilities, 
one still has uniformly bounded solutions if d < 1. This fact is easy to demonstrate by a 
simple comparison argument, Martin and Pierre [17]. However, when d > 1, no comparison 
argument seems to apply, and a completely different approach is necessary. 

Our method is to seek local L p estimates in space by studying certain localized nonlinear 
functionals of solutions. Similar functionals appeared first in [18], and later in [13]. Since 
our solutions are only bounded in maximum norm and have no spatial decay at infinities, 
the functionals in [18] and [13] are not directly applicable. As in Collet and Eckmann [12], 
and Collet [11], we introduce smooth cut-off functions and convert the global functionals 
of previous authors into local ones. The first kind of cut-off functions we employ are 
simply: if = <f(x) — (1 + \x — x \ 2 )~ n , where x is an arbitrary point in R n , and is used 
to translate the location of cut-off so as to achieve uniform L°° bounds in space. With 
such a if and the resulting local W estimates, we prove the existence of global classical 
solutions. However, the L°° norm of solutions grow exponentially in time. To improve 
the L°° estimates of solutions, we consider the second kind of cut-off functions which are 
time dependent solutions of the backward heat equation ip t + dA<f = 0. Using these time 
dependent cut-off functions, we are able to refine the L°° estimates to the order of loglog 
growth for any space dimension. Thus the possible growth of u 2 is practically extremely 
hard to observe even if it exists. On the other hand, it remains an interesting problem to 
prove or disprove the uniform L°° bounds on u 2 . Our main result is: 

Theorem 1.1 Consider the thermal diffusive combustion system (1.1) with nonnegative 
initial data (u\, it 2 ) |t=o — (^1,02) £ (C& U (-R n )) 2 ; an d d > 1. Then there exist unique global 
in time classical solutions (u u u 2 ) E C([0, +00); (C^ u (R n )) 2 ) n ^((0, +00); (C^ u (R n )) 2 ). 
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Moreover, let 11(01,02)1100 = max(| |ai| |oo, 11021100), then there exists positive constant C = 
C(n,d,m) such that: 

|K(t,x)|| L °°(i?™) < IKHoo, 

IMM)IU~(^) < C\\(a 1 ,a 2 )\\ 00 loglog(\\(a 1 ,a 2 )\\Zt + 2e). (1.2) 

Corollary 1.1 Under the same assumptions in the above theorem except that the nonlinear 
term uxu™ is replaced by either u\exp{ — ^} or u\u™exp{— E > 0, then there exists 
positive constant C = C(n,d, ||(ai, O2)||oo, m) such that: 

\\Ul{t,x)\\ L <x>( R n) < | |ai j |oo, 

\\u 2 (t,x)\\ L <x>( R n) < Cloglogit + 2e). (1.3) 

Remark 1.1 In case of power nonlinearities u\u™ , the system has the following scale 
invariance property: if ui = Ui(t,x), % = 1,2, are solutions, then so are Vi = Vi(t,x) = 
A~Mj(A 2 t, \x), for any A > 0. That is why the estimates are also in the scale invariace 
form in the theorem. In case of the Arrhenius reactions, we lose the scale invariance due 
to the exponential term exp{ — ^}, hence we do not know the explicit dependence of C on 
1 1 (oij 02)| |oo- The proofs are the same in both cases except for some technical details that 
we will point out later. 

Remark 1.2 // the initial data for u 2 , i.e. the function a 2 is strictly above any positive 
constant, then maximum principle shows that u 2 is above this constant forever, and so u\ 
decays to zero exponentially fast. By Theorem 1.1, u 2 is uniformly bounded for all time 
since uxu™ decays exponentially in time. 

The rest of the paper is organized as follows. In section 2, we use the first kind of 
time independent cut-off functions and local nonlinear functionals to prove the existence 
of global solutions. In section 3, we employ the time dependent cut-off functions and their 
properties to refine the L°° estimates of solutions and complete the proof of the main 
theorem. 

2 Global Existence of Classical Solutions 

In this section, we establish the global existence of the classical solutions of the thermal 
diffusive system: 
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«2,t = dA x u 2 + UiU™, (2.1) 

where x G -R™, t > 0, n > 1, m > 1, d > 1; the initial data (ai(x), a 2 (x)) are bounded 
uniformly continuous functions on R n , denoted by C® u (R n ). Local existence of nonnegative 
classical solutions on a maximal existence interval [0, T ) is standard, and we only need to 
derive estimates of solutions independent of T , so as to continue the classical solutions 
forever in time. 

Consider classical solutions u, t G C([0, T ); C$ u (R n )) n ^((0, T ); C b ° J, i = 1,2, for some 
T > 0. Let F = F(u 1 ,u 2 ) be a smooth function of Mj, such that F > 0, Fj > 0, 
> 0, i = 1, 2, here we abbreviate Fj = J^, similarly for the second derivatives. Let also 
y? = <f(t,x) be a smooth nonnegative function with exponential spatial decay at infinity. 
Writing / in place of J Rn ■ ■ ■ dx, we calculate using (2.1) and integration by parts: 

d t J ' <pF = J <p t F + J (pFiU^t + J (fF 2 u 2 ,t 

= J °ip t F + J <pF 1 Au 1 + dj <pF 2 Au 2 - J ip{F l - F 2 ) Ul u™ 
= J(p t F- J y?F l5l |VMi| 2 - (d + 1) j <pF 1>2 V Ul ■ Vu 2 - d J pF 2j2 \Vu 2 \ 2 
- J Fi Vy? • Vui - d J F 2 Vy? -Vu 2 - j ip{F x - F 2 )u 1 u 2 n . (2.2) 

In view of: 

- J F 1 Vy? • Vui - y F 2 Vy? • Vw 2 = y Ay?F, 

we get: 

d t J^F = J (<p t + dA v )F + (d - 1) y Fx Vv? • V«i 

- y ^[F^ilVuil 2 + (1 + d)F lt2 Vu 1 ■ Vu 2 + dF 2t2 \Vu 2 \ 2 } 

- y i P {F 1 - F 2 ) Ul uZ, (2.3) 

which is our basic identity. 

By maximum principle, ||ui||oo(0 — ll a i||oo < C < +oo; Ui > 0, with strict inequality for 
t > 0, any x G i? n . To apply (2.3), we require: 

Fi > 2F 2 , 

(d+l^^dFx,^, (2.4) 
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for any (1x1,1*2) £ [0, C] x Under the conditions (2.4), we have from (2.3): 
d t JvF< J (ip t + dA<p)F + (d - 1) y F x Vy? • Vixi 

- i 1 p[Fi,i|Vui| 2 + gLF 2i2 |Vw 2 | 2 ] 

- l -J V F lUl u™. (2.5) 
As a first application of (2.5), we choose: 

v = * {x) = (i + \x- Xo \*r 

F(u 1 ,u 2 ) = (A + u 1 + ul)e €u2 , (2.6) 

where x is an arbitrary point in R n so that we can translate the function ip to achieve 
uniform estimates of solutions in space; and A, e _1 suitably large to be determined. We 
verify conditions (2.4) as follows: 

F 1 = (l + 2 Ul )e eu \ 
F 2 = e(A + Ul + u\)e eu \ 



so: 



Also: 



Fi > 2F 2 , for u 2 ) e [0, C] x R 1 , if 2e{A + C + C 2 ) < 1. (2.7) 

F ltl = 2e eu \ 
F 2 , 2 = e 2 (A + u 1 + u 2 1 )e eu2 , 

F^ 2 = e{l + 2u 1 )e eu \ 



thus: 



{d+l?Fl 2 = (l + 2 Ul ) 2 (d+l) 2 < (l + 2C) 2 (rf + l) 2 ^ ] 
dF 11 F 2 2 2d(A + u 1 + u 2 1 ) ~ 2dA < ' 

if a > (i±jgm±i)!. (2 . 8) 

2a 

Combining (2.7) and (2.8), we see that for any given C and d, we first choose A according 
to (2.8) then e by (2.7). It follows from (2.5) for t E [0,T ): 

d t JpF< JdApF + (d- ^jFxVip- V«i - ^ J <fF hl \V Ul \ 2 . (2.9) 
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Now ip has the properties: 

\A<p\<K<p, \V<p\<K<p, (2.10) 
for some constant K > 0. We continue from (2.9): 

dtJtpF < JdK v F+ (d- l)K J F lV \V Ul \ - ^ J ^F^Vu^ 2 

< dKj V F+ l -(d-lfK 2 j V ^. (2.11) 

1 , 1 

Notice that: 

F 2 (\ _|_ Oil, "l 2 p e «2 

— < 2(A + Mi + ^)e e " 2 = 2F, since A > 1. (2.12) 



F 1A 2 
Finally we end up with: 

d t JpF< [dK + (d- 1) 2 K 2 ] J V F, 



or: 



J y?F < Te a \ (2.13) 

where 

a = dK + (d - 1) 2 K 2 , r < c(n)(A + C + C 2 )e^™, (2.14) 
where c(n) is a positive dimensional constant. 

We now prove that inequality (2.13) implies that for any unit cube Q and any finite p > 1: 

/ M p < A- 1 Te at 2 n e- p (p+ l) p+1 . (2.15) 
In fact, we have with any nonegative integer k: 

e»T > f vF > -4/,e- > > (2,6) 
by taking rr at the center of Q. Hence, 

/ \u 2 \ k < A~ l Ye ut Te- k k\, (2.17) 

JQ 

which implies (2.15) by interpolation. 
Now we look for L°° bound on w 2 - Let 

G T (z) = {A7idT)-^e-^, z G R n , 
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then 



u 2 (t, x) = G t * a 2 + f G t _ s * ( Ul i%)(s)ds. (2.18) 
Jo 

The first term on the right hand side of (2.18) is bounded by 11021100. Moreover, 

G t - s *( Ul vZ)(s,x) = (47rd(t-s))-3 J e-^(u 1 u 2 n )(s,y)dy. 

Let {Qj}, j — 0, 1, 2, • • •, be the tiling of R n by unit cubes Q/s such that x is at the center 
of Qq. We have: 

j e '^^(u 1 u 2 n )(s,y)dy = ]T f e ~^A ' (u 1 u 2 n )(s,y)dy. (2.19) 



For y G Qj, we have the inequality: 



dist(x ,Qj ) 2 



g 8d(t-s) < gup g 8d(i-s) — g 8d(t-s) 

y&Qj 

Also there exists positive dimensional constant c = c(n) such that if y e Qj, j 7^ 0, we 
have: 

c(n) dist 2 (x, Qj) > \x - y\ 2 . (2.20) 
Applying Holder's inequality with p > | and its conjugate g, we get: 

I e~&i(u 1 u 2 n )(s,y)dy < (/ e"wS d y )V* . ( / («?^)( a ,y)dy) 1 /i' 

< c(n)d n ^(t-sr^([ (u{uT)(s,y)dy) l lv 

< (t- s) n/29 fi(n,d,a i? t), (2.21) 

where c(n) is a dimensional constant and by (2.15): 

Q(n,d,a h t) = c(n)d n/2q \\a 1 \\ O0 [A~ 1 Te at 2 n e~ mp (mp+l) mp+1 ] 1/p 

= c(n,d,ai)e at/p e- m (mp + l) m+1/p , (2.22) 

here c(n,d,a,i) is a positive constant depending on n, d, and | |(ai, a 2 )| |<x>- 
We deduce from (2.18) (2.22) that: 

dist(x,Qj) p (x-y) 2 

G t _ s *«iu!?(s,2;) < (Audit- s)T n/2 5> sdjt^T) / e"53F : 5y( WlU ^)(s,j/)dy 

Qj jQ i 



< {4nd)' n/2 Q{n, d, <n, t)(t - s)- n/2p £ e _± W=¥ 

< {4nd)- n/2 Q{n,d,a t ,t){t-s)- n/2p {[ e"^)M cfy + 1) 

< (47rrf)-" /2 fi(n, d, a*, *)(* - s)- n/2p (c(n, - s) n/2 + 1) 

< c(n,d)fi(n,d,a J ,t)((t-s) n/29 + (t-s)" n/2p )). (2.23) 

So integrating (2.23) on s G [0,t] gives: 

n I -i i re 

| M*, a;) | |oo < |ja 2 j|oo + c(n, d)n(ra, d, Oi,f)(f a « + t 2 ^), (2.24) 

where p > n/2. Estimate (2.24) and the standard parabolic regularity theory then implies 
the global classical solution («i,u 2 )(*,a:) G (C([0, +oo); C° b ) n C 1 ((0, +oo); C° 5 )) 2 . 

3 Large Time Asymptotic Bounds of Solutions 

In this section, we improve the L°° estimates (2.24) from exponentially growing in time to 
the order of loglog growth and complete the proof of the main theorem. Let us choose as 
before: 

F = F(u u u 2 ) = (A + Ul + u\)e eu \ 
yet the function ip is now a solution of the backward heat equation: 

ip t + dAip = 0. 

Define: 

g(u) = A + u + u 2 , 

and 

g = u + u 2 . 

Let us consider t G [0, T), where T is a suitably large but fixed time. The function is 
explicit: 

ip = (p(t,T;x) = (4ird)-?(T-t)-?e~' i $=V. (3.1) 
With the above choice, inequality (2.5) gives: 

d t J<pF<(d-l) J F 1 (V<p-Vu 1 )-^ J <p[F 1:1 \V Ul \ 2 + dF 2:2 \Vu 2 \ 2 }-^ J <pF lUl u™. (3.2) 
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The first integral of the right hand side of (3.2) can be transformed using integration by 
parts as follows: 

J F^Vip ■ V«i) = J #'(«i)e eu2 (Vy9 • V«i) = J e eu2 Vy? • Vg( Ul ) 

= J e e " 2 Vy? ■ V£(«i) = - j Aipe eu2 g( Ul ) - e J e eua £(ui) ■ Vu 2 

= J1 + J2. (3.3) 

In view of (3.1), we see that: 

Ay, = -d"V, = -d-\4nd)-^(T - t)~^ - MT 'j 2 t)2+f ]e-^) 
n 1 __n!_ 

> P id(T-t) 

2d(47rd)"/ 2 (T-t)§ +1 



c(n, d) 
T - 

which implies that: 



> (3-4) 



Ji < fzf / <P9Me«». (3.5) 



|./ 2 | < / g 2 ^' V ^V 2 + i! /^ M | Vw2 | 2 e e -. (3.6) 



On the other hand, 



Combining (3.2) (3.6), we get (F 2)2 = e 2 #(ui)e e " 2 ): 



9 t |^F < (d - lfj^ J V ~g( Ul )e^ + (d - I) J 



c(», d) + - ✓(«,)*<]. (37) 



Notice that: 



„ 1 1 -X M 2 

^ (47rd)™/ 2 (T-t)f 2d(T-t) 

= -c(n,d)— -<p, 

|Vy?| 2 c(n,d) |:r| 2 
V? 2 ~ (T-t) 2 ' 



(3.8) 
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which along with (3.7) yields: 

a t j, F < C( „, d) /^ [ HL±^ + ( Jl + _;^f_ t)2 -(i + 2Ml) ^l 

< c(n, d, | ^e eu2 u x {^— t + ^ - <) 

/■ | |x| 2 "U m 

< c(n,d, ai ) rn y?e e " 2 m(— — + ' , 2 - -f) 
7 {xe/? n| + <_L_ + j^ } T-t (T-£) 2 2 



2 -T-t T (T-t)2 



c(n, d, Oj 



Jipe eu2 u lU ^. (3.9) 



Remark 3.1 //we /iave £/ie Arrhenius type reactions g ar = u™exp{ — ^} , which is mono- 
tonely increasing in w 2; we simply need to replace the explicit inverse (gar)" 1 by an implicit 
one g~} . Everything else is the same. 

If T — £ > 1, then requiring e < ^ when m = 1, we get: 

c(n,d,ai) /\ . « m 2 . Id 2 . 2 er=J-+-L^)™ 



5, 



/• „ c(n, d, aj /\ . n w ir , 2e(-*—+ w , 

JVF < -7f±f- J(T-t)~e «(^i) ( i + _LL_ )e (T-t+pH^ 



c(n,d,ai) f r „„ 

— - y v?e 2 -ui«™ 



c(n,d,ai) /\ , „ _ M 2 . bl 2 , , f / u M 2 \ 5 
^ V_£ - J(T-t)-*e i^ (1 + _L_L_ )e 26(i+— ). 



c(n,rf, aj) c(n,rf, aj) /■ eu2 m 

" T-£ 2 MlM2 ' 

which implies that: 

°^ y* y v?e eu2 Mi< + J^F < J (p(0, T; x)F{a u a 2 )dx + c(n, d)Zo$T 

= ||-^(oi, 02) | |oo + c(n, d)logT, (3.10) 

for £ G [0, T — 1). It follows from (3.10) that: 

[ T 1 dt [ ipe^u^dx < c(n, d, a;)(l + logT) (3.11) 

JO JR n 
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By choosing (p = ip(t, T;x — x ), for any x G R n , we also arrive at (3.11) and so: 

[ T 1 c%*(e eu2 wi<)(t) < c(n,d,ai)(l + logT). (3.12) 

JO 

Let v = u 2 , k > 2, then 

v t = ku^~ 1 u 2 , t , 
= ku 2 ' 1 Vu 2 , 
Av = fcu^Aua + fc(ib - 1)m(T 2 |Vw 2 | 2 , 
and so v satisfies the equation: 

v t = dAv - dk{k - 1)m^ 2 |Vm 2 | 2 + kuiu 2 n+k ~ 1 , (3.13) 

which implies: 

v t <Gt*v + k f dsGts + ium^™- 1 )^^). (3.14) 
J o 

Letting t = T in (3.14) yields: 

< ||a 2 ||* +Jfe f ^ dsG T _ s * (u 1 u k 2 +m ~ 1 )(s) + k f dsG T _ s * (uju^-^fs), 
jo " ir-i 

which shows by (3.12): 

u(T,x)< | |a 2 1 lSo + c ( n > d, a i )e~ k k\(l + Zo<?T) + k [ dsG T ^ s * (u l u k 2 +m - 1 )(s). (3.15) 

Jt-i 

Now it suffices to consider: 

J V (t,T;x) Ul u k 2 +m -\x), 

where t G [T - 1,T]. 

Let {Qj}^ be a tiling of R n with unit cubes, and located at the center of Q . Then: 
J tp(t, T; x)uiu k , +m ~ 1 (x) < c(n,d)£ J Q (y e'^mu^dx 

< c(n,d)^— — -^e 8d ( T -') ^ e'^^MiM^™" 1 ^. (3.16) 

Following the argument in (2.19) and (2.21), we get: 

f e-^^u^^dx <c(n,d)(T -t)%([ ul {k+m ~ 1] {t,y)dy)* . (3.17) 

JQj JQ 3 
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On the other hand, adjusting T to T + 1 in (3.10), we have for t e [T - 1, T]: 

J ip(t,T + l;x)F(u 1 ,u 2 )dx < c(n,d, Oj)(l + log(T + 1)), 

or 

J e -«(T+i- t) e ^dx < c{n, d, ai )(l + logT), 

or 

J e-^fe €u2 dx < c(n, d, a;)(l + logT). 
Using the spatially translated <p, we find: 

/ e eu2 dx < c(n, d, Oj)(l + logT), V j, 

or 

/ w^dx < 7!c(n, d, ai)e - T(l + %T), V 7, j, e Z + , i e [T - 1, T). (3.18) 
By (3.18), (3.17) gives by Holder's inequality: 

[ e~^^u x u\ +m - 1 dx<c{n,d,a i ){T-t)^{t uf +1 dy)*M+v, 

JQj JQj 

n 3 3 3 

< c(n,d,ai)(T -t)^(([f3] + l)!)5n3J+U e "(l + logT)^m+ri y (3.19) 
where (3 = p(k + m — 1), stands for the integral part of (3. 

Proof of Theorem 1.1: Combining (3.16) and (3.19) and again following the argument 
in (2.19) (2.24) shows for T > e: 

J V (s,T;x)u l u k 2 +m ~ 1 {s,x)dx 

n n 3 3 3 

< c{n,d, ai ){(T - s)~« + (T - s)"^)(([/3] + l)\)im+Ve-*(logT)*M+V, (3.20) 
where p > max(l, |). Integrating (3.20) from T — 1 to T shows via (3.15) that: 

v(T,x) < ||a 2 ||So + c(n,d, ai )e- k (logT)k\ + c(n,d,a i )e~?((p(k + m))\) 1/p (logT) 1/p , 

or: 

u 2 (T, x) < I \a 2 \ loo + c(n, d, a t )(k\logT + [logflf/ p {{p{k + m))!) 1 ^) 1 ^, 
by Sterling's formula ( + to))! < c(n,m)e pkl ° 9k ): 

< C(n, d, Oi, m )( e Mogk e lo g logT + e loglogT e klogky/k 
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< c(n, d, a h m)e lJ ^ si:+lo9k . (3.21) 
Minimizing the exponent with respect to k shows that we should choose: 

k = [loglogT] + 1, (3.22) 

which implies from (3.21), and (2.24) that: 

U2(T, x) < c(n, d, di, m)loglog{T + 2e) (3.23) 

for all T > 0, where constant c > depends on n, d, ||(ai, calico, and m. In case of power 
nonlinearities uiu™, we first consider initial data such that ||(ai, Q2)||oo < 1 and so drop 
the dependence on in the bound (3.23). We verify by direct substitution that if Ui(t,x), 
% = 1, 2, are solutions, then for any A > 0, \™Ui(\ 2 t, Xx) are also solutions. Now choose A 
such that 

A"£ = IKai.aaJlloo. (3.24) 

— 9 

It follows that the L°° norm of the initial data of the solutions A m Wj(A t, Xx) is equal to 
one. By (3.23), we have: 

A™ | |(-ui, -u 2 )(A t, Ax) | |oo < c(n,d,m)loglog(t + 2e), (3.25) 

which implies: 

| \(ui, u 2 )(X 2 t, x)\ |oo < c(n, d, m)X~™loglog(t + 2e). 
Rescaling time t, we obtain: 

1 1 (tii, u 2 )(t, x)\ |oo < c(n, d, m)X~™loglog(X~ 2 t + 2e), 

which is just: 

1 1 (u u u 2 )(t, x) | |oo < c(n, d, m) 1 1 (ai, a 2 ) | U%MI I a 2 ) 1 1™* + 2e), (3.26) 
by recalling (3.24). The proof of the main theorem is complete. 
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